Current in a three-dimensional periodic tube with unbiased forces 
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- - - Transport of a Brownian particle moving along the axis of a three-dimensional asymmetric periodic 

' tube is investigated in the presence of asymmetric unbiased forces. The reduction of the coordinates 

' may involve not only the appearance of entropic barrier but also the effective diffusion coefficient. 

■ ft is found that in the presence of entropic barrier, the asymmetry of the tube shape and the 

asymmetry of the unbiased forces are the two ways of inducing a net current. The current is a 
I peaked function of temperature which indicates that the thermal noise may facilitate the transport 

' even in the presence of entropic barrier. There exists an optimized radius at the bottleneck at which 

the current takes its maximum value. Competition between the two opposite driving factors may 
induce current reversal. 

Q ■ PACS numbers: 05. 60. Cd, 05. 40. Jc, 02. 50. Ey 
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Transport phenomena play a crucial role in many processes from physical, biological to social systems. There has 

Ci • been an increasing interest in transport properties of nonlinear systems vi^hich can extract usable vifork from unbiased 

^ nonequilibrium fluctuations [H, 13) S S iB] ■ In these systems, directed Brownian motion of particles is generated by 

' nonequilibrium noise in the absence of any net macroscopic forces and potential gradients ,4]. In all these studies, the 

^ , potential is taken to be asymmetric in space. It has also been shown that a unidirectional current can also appear 

Q ' for spatially symmetric potentials if there exits an external random force either asymmetric [2] or spatially-dependent 

The most studies have been referred to the consideration of the energy barrier. The nature of the barrier depends 
, on which thermodynamic potential varies when passing from one well to the other, and their presence plays an 
^ important role in the dynamics of the solid state physics system. However, in some cases, such as soft condensed 
matter and biological systems, the entropy barriers should be considered. The entropy barriers may appear when 
coarsening the description of a complex system for simplifying its dynamics. Reguera and co-workers Q use the 
mesoscopic nonequilibrium thermodynamics theory to derive the general kinetic equation of a system and analyze 
in detail the case of diffusion in a domain of irregular geometry in which the presence of the boundaries induces an 
entropy barrier when approaching the dynamics by a coarsening of the description. In their recent work 7] , they study 
the current and diffusion of a Brownian particle in a symmetric channel with a biased external force. They found 
, that the temperature dictates the strength of the entorpic potential, and thus an increasing of temperature leads to 
a reduction of the current. The previous works on entropic barriers is limited to biased external force. However, in 
some cases, the net current occurs in the absence of any net macroscopic forces or in the presence of unbiased forces. 
For example. Molecular motors move along the microtubule without any net macroscopic force. The present work is 
^ , extend the study of entropic barriers to case of unbiased forces that indicates zero force at macroscopic scale. Our 
- - ' emphasis is on finding conditions of obtaining a net current in the presence of entropic barriers. The asymmetry of 
the tube shape and the asymmetry of the unbiased fluctuations are the two driving factors for obtaining a net current. 
When the two driving factors compete each other, the current may reverse its direction. 

CURRENT IN A THREE-DIMENSIONAL PERIODIC TUBE 

We consider a Brownian particle moving in a asymmetric periodic tube [Fig. 1] in the presence of an asymmetric 
unbiased fluctuations. Its over-damped dynamics is described by Q 
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FIG. 1: Schematic diagram of a tube with periodicity L. The shape is described by the radius of the tube u){x) — a[sin(2^) + 
■J sin(i^)] + b. A is the asymmetric parameter of the tube shape. 



where r is the 3D coordinate, rj the friction coefficient of the particle, fcs the Boltzmann constant, T the temperature, 
and ^ (t) is Gaussian whit noise with zero mean and correlation function: < [t ) 25ij5{t—t ) for i, j = a;, y, z. 

< ... > denotes an ensemble average over the distribution of ^ (t). S{t) is the Dirac delta function. The reflecting 
boundary conditions ensure the confinement of the dynamics with the tube. Fx{t) is an asymmetric unbiased external 
force along x direction. F{t) is its scalar quantity and satisfies [1, Q 

F{t + T) = F{t), f F{t)dt^O, (2) 
Jo 



^ ' \~Fo, nT+^T{l-e) <t<{n+l)T, ^•'> 

where r is the period of the unbiased force , Fq its magnitude and e the temporal asymmetric parameter with 
-1 < e < 1. 

The shape of the tube is described by its radius 

, , , . .2nx. A . Anx... , 
^(x) = a[sm(— ) + - sm(— )] + b, (4) 

where a is the parameter that controls the slope of the tube, A the asymmetry parameter of the tube shape. The 



radius at the bottleneck is r;, = 6 — ay 1 + 

The movement equation of a Brownian particle moving along the axis of the 3D tube can be described by the 
Fick-Jacobs equation 0, 0, [13) Hi which is derived from the 3D (or 2D) Smoluchowski equation after elimination 
of y and z coordinates by assuming equilibrium in the orthogonal directions. The reduction of the coordinates may 
involve not only the appearance of entropic barrier, but also the effective diffusion coefficient. When \uj {x)\ << 1, 
the effective diffusion coefficient reads 

[1 + W (X)'']'^ 

where Dq — ksT/rj and a = 1/3, 1/2 for 2D and 3D, respectively. 

Consider the effective diffusion coefficient and the entropic barrier, the dynamics of a Brownian particle moving 
along the axis of the 3D tube can be described by @, 0] 
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where we defines a free energy A{x,t) := E — TS = —F(t)x — TkB^'n.h{x), here E = —F{t)x is the energy, S = 
kB^n.h{x) the entropy, h{x) the dimensionless width 2uj{x)/L in 2D, and the dimensionless transverse cross section 
7r[a;(x)/L]^ of the tube in 3D. j{x,t) is the probabihty current density. P{x,t) is the probabihty density for the 
particle at position x and at time t. It satisfies the normalization condition P{x,t)dx — 1 and the periodicity 
condition P{x, t) = P{x + L,t). 

If F{t) changes very slowly with respect to t, namely, its period is longer than any other time scale of the system, 
there exists a quasi-steady state. In this case, by following the method in [l|, 0, 0, |j, 0, [^, 0, , we can obtain the 
current 



Jim) 



kBT[l — exp(- 



lo ^(2^)exp( 



(7) 



The average current is 
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with 



Ji = (1 - e)j{\^Fo), j2 = (1 + e)j{-Fo). 
1 — e 



(9) 



RESULTS AND DISCUSSIONS 

Because the results from 2D and 3D are very similar, for the convenience of physical discussion, we now mainly 
investigated the current in 3D with fcs = 1 and rj = 1. 
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FIG. 2: Current J versus temperature T for different asymmetric parameters A of the tube shape at a = 1/2tt, b — 1.5/27r, 



L = 2tt, a = 1/2, Fo = 0.5 and e = 0.0. 



The current J as a function of temperature T for the case e = 0.0 is presented in Fig. 2 for different values 
of asymmetric parameter of the tube shape. The curve is observed to be bell shaped, which shows the feature of 
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resonance. When T ^ 0, the particle cannot reach 3-D area and the effect of cntropic barrier disappears and there is 
no current. When T ^ oo, the effect of the unbiased forces disappears and the current goes to zero, also. There is 
an optimized value of T at which the current J takes its maximum value, which indicates that the thermal noise may 
facilitate the particle transport even in the presence of cntropic barrier. The current is negative for A < 0, zero at 
A = 0, and positive for A > 0. Therefore, the asymmetry of the tube shape is a way of inducing a net current. 




FIG. 3: Current J versus the radius ri at the bottleneck at o = l/27r, a = 1/2, L = 2-k, Fq = 0.5, T = 0.5, A = 1.0 and 
e = 0.0. 

Figure 3 shows the current J as a function of the radius at the bottleneck n, for the case £ = 0.0 and A = 1.0. If 
the bottleneck has zero the particle can not pass through the bottleneck, so the current should be zero. When the 
bottleneck has infinite radius, the effect of tube shape disappears and the current tends to zero, also. Therefore, the 
current J is a peaked function of the radius rt, at the bottleneck. 

Figure 4 shows the current J versus the amplitude Fq of the unbiased forces of for the case e = 0.0 and A = 1.0. 
When Fq ^ 0, only the effect of the cntropic barrier exists, so the current tends to zero. The current J saturates to 
a certain value in large amplitude Fq limit. There exists an optimized value of Fq at which the current T takes its 
maximum value. 

The current J as a function of temperature T for the case A = 0.0 is presented in Fig. 5 for different values of 
asymmetric parameter e of unbiased forces. The curve is a bell shaped function of temperature which is the same as 
that in Fig. 2. The asymmetry of the unbiased forces can induce a net current even when the shape of the tube is 
symmetric. For the asymmetry of the tube shape, the asymmetry of the external force is another way of obtaining a 
net current. Similarly, J is negative for e < 0, zero at £ = and positive for £ > 0. 

In Fig. 6, we plot the current J as a function of temperature T for different combinations of A and s. When the 
asymmetric parameter A of the tube shape is positive, the current may reversal its direction on increasing temperature 
for negative e. It is obvious that the current reversal may occur when a positive driving factor competes with a negative 
one. Therefore, there may exist current reversal for Ae < (A = 1.0, s = —0.1; A = 1.0, e = —0.15; A = 1.0, 
£ = —0.2). However, A£ < is not a sufficient condition for current reversal. For example, the current is always 
positive for A = 1.0, £ = —0.05. 




FIG. 5: Current J versus temperature T for different asymmetric parameters e of the external force at o = l/27r, b = 1.5/27r, 
L = 2-K,a= 1/2, Fo = 0.5 and A = 0.0. 

CONCLUDING REMARKS 



In this paper, we study the transport of a Brownian particle moving in a 3D periodic tube in the presence of the 
unbiased forces. The presence of the boundaries induces an entropic barrier when approaching the exact dynamics by 
coarsening of the description. Both the asymmetry of the tube shape and the asymmetry of the unbiased forces are 
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FIG. 6: Current J versus temperature T for different combinations of e and A at a = l/27r, b — 1.5/27r, L = 2tt, a — 1/2 and 



Fo = 0.5. 



the two driving factors for obtaining a net current. When the two driving factors compete each other, the current may 
reverse its direction upon increasing temperature. When the bottleneck has zero radius and infinite radius then the 
current tends to zero. An optimized bottleneck radius leads to a maximum current. There exists an optimize value 
of temperature at which the current takes its maximum value, which indicates that the thermal noise may facilitate 
the current. The results we have presented have a wide application in many processes [6|, such as molecular motors 
movement through the microtubule in the absence of any net macroscopic forces jl3| . ions transport through ion 
channels [3], motion of polymers subjected to rigid constraints 15], drug release [lit and polymer crystallization [l3]. 
In these systems, a directed-transport with entropic barriers can be obtained in the absence of any net macroscopic 
forces or in the presence of the unbiased forces. 
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